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Abstract 

The almost frictionless transport of the very-high temperature amorphous matter which 
resembles the color glass condensate (possibly having much of their origin in the RHIC heavy 
ion collisions) in a confined annular tube with transversely corrugations is investigated by us- 
ing the verified transition-rate model and boundary perturbation method. We found that for 
certain activation volume and energy there exist possible frictionless states which might be 
associated with the perfect fluid formation during the early expansion stage in RHIC Au+Au 
collisions. We also address the possible similar scenario in LHC Pb+Pb collisions considering 
the possible perfect fluid formation in ultra-high temperature transport of amorphous matter. 

KEYWORDS : Phenomenological Models, Statistical Methods, Quantum Dissipative Sys- 
tems, Thermal Field Theory, Boundary perturbation. 

1 Introduction 

The success of ideal hydrodynamics [1-7], in explaining bulk of the data in Au+Au collisions 
at RHIC (cf., e.g., [8-9]), has led to a paradigm that in Au+Au collisions a nearly perfect 
fluid is created. There is no a priori reason why ideal hydrodynamics works so well at RHIC. 
The relation between the inverse Reynolds number (-R -1 ) and the attenuation length (L s = 
(4/3)ry/(e + p)) becomes R^ 1 ~ L s /r in the Bjorken flow case (rj is the shear viscosity, e is the 
energy density, p is the pressure, r(~ 0.6 fm/c) is the expansion rate) [10-11]. FT 1 is found to 
be very small from the blast-wave fitting with the above correction term (cf. [10]). This suggests 
that the hadronic fluid in heavy ion collisions is nearly a perfect one. 

The perfect fluid can also be understood from the form of the stress-energy tensor. In the 
approximation where the Lagrangian is C = P(X), with P (X) being the kinetic function, the 
stress-energy tensor has the form 

oc -P(X) 9llu + 2P'(X) Ufl u v , (1) 

where = d^(p = u > <ft 1S a (real) scalar. This has the form of a stress-energy tensor for a 
perfect fluid with 4- velocity [12]. The fact that is a gradient means that the flow is 
irrotational : V x u = 0. Note that the defining property of a perfect fluid (fluid with zero 
viscosity [13]) is that at each point there is a frame in which the stress-energy tensor has the 
form T^ u = diag(— p,p,p,p) and in ideal hydrodynamics, the mean free path of particles is 
assumed to be zero. 
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The interesting (relevant to our borrowed approach here) issue is the possible existence of the 
color glass condensate (CGC) (cf. e.g., [13,14,7,8]). Au+Au collisions are then collisions of two 
sheets of colored glass, with the produced quarks and gluons materializing at a time given by the 
inverse of the saturation momentum. The proposed color glass condensate was thought to be a 
possible precursor state to the quark-gluon plasma (QGP). To be specific, at very high energies, 
multi-particle production in QCD is generated by small x partons in the nuclear wavefunctions 
(x is the fraction of the longitudinal momentum carried by the parton). These partons have 
properties best described as a color glass condensate (CGC) [5]. When two sheets of CGC collide 
in a high-energy heavy-ion collision, these partons are released and create energy densities an 
order of magnitude above the energy density required for the crossover from hadronic to partonic 
degrees of freedom. This matter, at early times after a heavy-ion collision, is a coherent classical 
field, which expands, decays into nearly on shell partons and may eventually thermalize to form 
a quark-gluon plasma. Because it is formed by melting the frozen CGC degrees of freedom 
and is the non-equilibrium matter preceding the QGP, this matter is called the glasma [5,13] (a 
non-equilibrium gluonic state between the collision moment and equilibrated QGP [6]). 
As remarked in [13], the color glass condensate is the matter associated with the wee partons of a 
high energy hadronic wavefunction. This matter has very high energy density. It has properties 
similar to Bose-Einstein condensates (there is a characteristic momentum scale, termed the 
saturation scale, below which the gluon density saturates. This effect sets in when x becomes 
small and the associated gluon wave length increases to nuclear dimensions. In such a regime 
gluons may interact and form a coherent state reminiscent of a Bose-Einstein condensate) and 
to spin glasses. It was motivated based on observations in deep inelastic scattering. The idea is 
that as one goes to higher and higher energies, there are more and more gluons in the hadron 
wavefunction. They have to squeeze together, and highly occupy phase space, making a highly 
coherent high density system of gluons. This matter controls the high energy limit of hadronic 
scattering and provides the initial conditions for the matter made in such collisions. Thus 
the almost perfect fluid phenomena (or frictionless transport) observed during the early stage 
(expansion) for Au+Au collisions in RHIC could be dominated by small x (gluons) (so that the 
geometrical scaling occurs) [13,14]. Note that the gluon population at low x is not an incoherent 
superposition of nucleon structure functions but is limited with increasing Au by nonlinear 
gluon-gluon fusion resulting from the overlap of gluons from several nucleons in the plane of the 
nucleus transverse to the collision axis. 

In fact due to the nearly perfect fluid (very low viscosity or flow resistance being almost zero) 
characteristics associated to a possible new state of matter created in Au+Au collisions at 
RHIC some researchers argued that this property (the (kinematic) shear viscosity of this nearly 
perfect fluid has been determined and found to be at least a factor of 4 smaller than that of 
the superfluid 4 He [2]) could be called high temperature superfluidity because the fluid of quarks 
in the super-high temperature range flows better than the superfluid 4 He in the extremely low 
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temperature region (cf. [2,4]). With above facts we can have in mind the possible link between 
the (color) glassy condensate (with the role of gluons) and their almost frictionless transport 
(of the fluid or liquid). This latter situation is similar to the proposed superglass reported quite 
recently [15,16]. 

Meanwhile we know that at low temperatures gases condense into the liquid (or solid) state. 
In a liquid transport is no longer governed by the motion of individual molecules or composite 
(condensed) particles. As noticed in [1], Eyring proposed that momentum transport of condensed 
composite particles or molecules (liquid-like) is due to processes that involve the motion of 
vacancies [17-19]. These processes can be viewed as thermally activated transitions in which 
a molecule or a cluster moves from one local energy minimum to another. Eyring's approach 
involves the Planck constant h and the appearance of h is related to Eyring's assumption that 
the collision time of the molecules is h/^ksT), the shortest timescale in a liquid (fee is the 
Boltzmann constant, T is the temperature). One important observation is : Eyring's approach 
has been successfully applied to the study of transport of glassy matter at a wide range of 
temperatures [17-19]. Possible frictionless transport of glassy matter under specific environment 
were also reported recently [16,20]. 

Above mentioned analogy or observed facts is the primary motivation for our present study. 
Our concern is associated with the possible superfluidity (or perfect fluid) formation after shear- 
thinning (i.e., the viscosity diminishes with increasing shear rate.). While superflow (or perfect 
fluid flow) in a state of matter possessing a shear modulus might initially seem untenable, 
experimental claims for precisely this phenomenon in solid 4 He now abound (cf. [15,16,20]. 
In this paper, to demonstrate the nearly frictionless transport of CGC and identify the possible 
transition (or critical) temperature, we shall adopt the verified transition-rate-state model [17- 
19] to study the transport of CGC (presumed to be amorphous) within a corrugated annular tube 
(shell-like) . The possible nearly frictionless states due to strong shear-thinning will be relevant 
to the perfect fluid formation at the early stage for Au+Au collidions in RHIC as mentioned 
above. To obtain the law of shear-thinning matter for explaining the too rapid annealing at the 
earliest time, because the relaxation at the beginning was steeper than could be explained by 
the bimolecular law, a hyperbolic sine law between the shear (strain) rate : 7 and shear stress : 
r has been proposed and the close agreement with experimental data was obtained [17-19]. This 
model has sound physical foundation from the thermal activation process [1,17-19] (a kind of 
quantum tunneling which relates to the matter rearranging by surmounting a potential energy 
barrier was discussed therein). With this model we can associate the (shear-thinning) fluid 
with the momentum transfer between neighboring clusters on the microscopic scale and reveals 
the possible microscopic interaction in the relaxation of flow with dissipation (the momentum 
transfer depends on the activation (shear) volume : V* = Vh which is associated with the center 
distance between atoms and is equal to /c^T/to (tq a constant with the dimension of stress). 
To consider the more realistic but complicated confined boundary conditions in the interfaces of 
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the annular tube (shell- like), however, we will use the boundary perturbation technique [16,20] to 
handle the presumed wavy-roughness along the interfaces of the confined annual tube. To obtain 
the analytical and approximate solutions, here, the roughness is only introduced in the radial 
or transverse direction. The relevant boundary conditions along the wavy-rough surfaces will 
be prescribed below. We shall describe our approach after this section : Introduction with the 
focus upon the transition-rate approach and boundary perturbation method. The approximate 
expression of the transport is then demonstrated at the end. Finally, we will illustrate our results 
into three figures and give discussions therein. 

2 Theoretical Formulations 

Researchers have been interested in the question of how (complex) matter responds to an ex- 
ternal mechanical load. External loads cause transport, in Newtonian or various types of non- 
Newtonian ways. Amorphous matter, composed of polymers, metals, or ceramics, can deform 
under mechanical loads, and the nature of the response to loads often dictates the choice of 
matter in various applications. The nature of all of these responses depends on both the tem- 
perature and loading rate [1]. 

To the best knowledge of the author, the simplest model that makes a prediction for the rate 
and temperature dependence of shear yielding is the rate-state model of stress-biased thermal 
activation [17-19]. Structural rearrangement is associated with a single energy barrier E that is 
lowered or raised linearly by an applied stress a : R± = vq exp[— E/iksT)] exp[±aV* / (hsT)], 
where vq is an attempt frequency and V* is a constant called the 'activation volume'. In amor- 
phous matter, the transition rates are negligible at zero stress. Thus, at finite stress one needs 
to consider only the rate i?+ of transitions in the direction aided by stress. 
The linear dependence will always correctly describe small changes in the barrier height, since 
it is simply the first term in the Taylor expansion of the barrier height as a function of load. 
It is thus appropriate when the barrier height changes only slightly before the system escapes 
the local energy minimum. This situation occurs at higher temperatures; for example, New- 
tonian transport is obtained in the rate-state model in the limit where the system experiences 
only small changes in the barrier height before thermally escaping the energy minimum. As the 
temperature decreases, larger changes in the barrier height occur before the system escapes the 
energy minimum (giving rise to, for example, non-Newtonian transport). In this regime, the 
linear dependence is not necessarily appropriate, and can lead to inaccurate modeling. To be 
precise, at low shear rates (7 < j c ), the system behaves as a power law shear-thinning material 
while, at high shear rates, the stress varies affinely with the shear rate. These two regimes 
correspond to two stable branches of stationary states, for which data obtained by imposing 
either a or 7 exactly superpose. 

We shall consider a steady transport of the amorphous matter (CGC) in a wavy-rough annu- 
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lar tube of 77 (mean-averaged inner radius) with the inner interface being a fixed wavy-rough 
surface : r = 77 + esin(/c# + (3) and r<i (mean-averaged outer radius) with the outer interface 
being a fixed wavy-rough surface : r = ri + esin(fc0), where e is the amplitude of the (wavy) 
roughness, (3 is the phase shift between two walls, and the roughness wave number : k = 2n/L 
(L is the wavelength of the surface modulation in transverse direction). 

Firstly, this amorphous matter (composed of composite condensed particles, say, quarks and 
gluons) can be expressed as [16,18,20] 7 = 70 sinh(r/ro), where 7 is the shear rate, r is the shear 
stress, ro = 2ksT/Vh, and 7o(= CkksT exp(—AE/kBT)/h) is with the dimension of the shear 
rate; here = 2Vf l /V m is a constant relating rate of strain to the jump frequency (V^ = A2A3A, 
V m = A2A3A1, A2A3 is the cross-section of the transport unit on which the shear stress acts, A 
is the distance jumped on each relaxation, Ai is the perpendicular distance between two neigh- 
boring layers of particles sliding past each other), accounting for the interchain co-operation 
required, h is the Planck constant, AE is the activation energy In fact, the force balance gives 
the shear stress at a radius r as r = —(rSQ)/2 [16,20]. 5Q is the net effective forcing along the 
transport (or tube-axis : z-axis) direction (considering dz element). 

Introducing the forcing parameter $ = — (r2/2ro)5Q then we have 7 = 70 sinh($r/r2). As 
7 = —du/dr (u is the velocity of the transport in the longitudinal (z-)direction of the annular 
tube), after integration, we obtain 

u = u s + ^[cosh*-cosh(— )], (2) 

here, u s (= u s n p ) is the velocity over the (inner or outer) surface of the annular (cosmic) string, 
which is determined by the boundary condition. We noticed that a general boundary condition 
for transport over a solid surface [16,20] was 

Su = L° 7 (l - |)- 1/2 , (3) 

where 5u is the velocity jump over the solid surface, L° s is a constant slip length, j c is the critical 
shear rate at which the slip length diverges. Note that the slip (velocity) boundary condition 
above (related to the slip length) is closely linked to the mean free path of the particles together 
with a geometry-dependent factor (in low temperature regime it is the quantum-mechanical 
scattering of Bogoliubov quasiparticles which is responsible for the loss of transverse momentum 
transfer to the confined interfaces [21]). The value of j c is a function of the corrugation of 
interfacial energy. 

With the slip boundary condition [16,20], we can derive the velocity fields and transport rates 
along the wavy-rough annular (cosmic) string below using the verified boundary perturbation 
technique [16,22] and dimensionless analysis. We firstly select LP S to be the characteristic length 
scale and set r' = r/L®, R\ = ri/L®, R2 = ^/L^, e' = e/L s . After this, for simplicity, we drop 
all the primes. It means, now, r, R±, R2 and e become dimensionless ($ and 7 also follow). The 
wavy interfaces are prescribed as r = R2 + e sin(A;^) and r = R\+e s'm(k6 + (3) and the presumed 
steady transport is along the z-direction (annulus-axis direction). 
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2.1 Boundary Perturbation 

Along the outer interface (the same treatment below could also be applied to the inner interface), 
we have 7 = (du)/(dn)\ on interfaces- Here, n means the normal. Let u be expanded in e : 

u = u + eui + e 2 n 2 H , 

and on the boundary, we expand u(ro + edr, 9(= 60)) into 

dv 2 

u(r,6)\( ro+edr>eo ) = u(r Q ,0) + e[dru r (r ,9)} + e 2 [—u rr {r ,9)} H = 

A /R 2 

{usiip + — [cosh$ - cosh( — )]}| on interfaces, r = Ri, R 2 ; (4) 

where 

u slip\on interfaces ~~ on interfaces' 

(5) 

Now, on the outer interface (cf., e.g., [16,22]) 

du _ V(r — R 2 — esm(k9)) M ok 2 2/ ,„s n _i r 
7 = Tn = VU ■ Mr-R 2 -e S in(k9))\ = [1 + 6 7* C ° S W1 

A; 

cos(/c6»)u e | (jR2+edrie) ] = u 0r |fi 2 +e[ui r | jR2 +u 0rr |/? 2 sin(fe^)- 
1 A: 2 

— uo e \n 2 cos(k9)] + e 2 [--— cos 2 (k9)u 0r \ R2 + u 2t \r 2 +u lrr \ R2 sm(k9)+ 

IT Zi T* 

1 k 

~ u o rrr \R 2 sm 2 (k9) - —cos(k9)(u le \ R2 + u 0er \R 2 sin(fc0))] + 0(e 3 ). (6) 

Considering L° s ~ i?i , i? 2 3> e case, we also presume sinh<E> <C j c /io- With equations (2) and 
(6), using the definition of 7, we can derive the velocity field (u) up to the second order : 

u (r,9) = -( J R 2 7o/^){cosh($r/i? 2 ) -cosh$[l + e 2 $ 2 sin 2 (M)/(2 J R|)] + 

e$sinh$ sm(k9)/R 2 } + u sHp \ r=R2+eshl ( ke y 

The key point is to firstly obtain the slip velocity along the boundaries or surfaces. After lengthy 
mathematical manipulations, we obtain the velocity fields (up to the second order) and then we 
can integrate them with respect to the cross-section to get the transport (volume flow) rate (Q, 
also up to the second order here) : 

rS P rR 2 +esin(k$) 
'Rl+esin(fc6»+/3) 

In fact, the approximate (up to the second order) net transport (volume flow) rate reads : 



fO p rR 2 +esm(k9) 

Q = / / u(r, 9)rdrd9 = Q + e Q po + e 2 Q P2 . 

JO JR 1 +esm(ke+l3) 



Q = nj {L° s (R 2 - i? 2 )sinh<D (1 - ^?)" 1 / 2 + ^[(R 2 - i? 2 )coshd> - |(i? 2 sinh<I>- 

7c/7o * * 
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i?!# 2 S inh($|i)) + ^|(cosh$ - cosh($|i))]} + e 2 {|n s%30 (i?| - i? 2 ) + 
L^7 sinh$(l + ^^)(-fc 2 + $ 2 )[1 - A 2 ] + ^7o[i?i shm^*) _ R 2 sinh <I>]- 

4 7c/70 -^2 ^ ^2 

-70 — [cosh<I> - cosh($ — )] + -70$ cosh <J>[i? 2 - 77-] + 
Z V XI2 4 XI2 

7T7o{[sinh$ + L°cosh<I>(l + Smh ^ )](i? 2 - i? 2 cos/3)} + -70^ cosh $+ 

7c/7o 2 5> 

r n 7r ^2- cosh$ , ^ /-n 

L °I $ ToTaT 1 "^) ]}cosh$. (7) 

4 7c/7o -n-2 



Here, 

= 2#Vo[sinh$(l - ^i)-V2]. (8) 

7c/ 7o 

3 Results and Discussions 

We firstly check the roughness effect (or combination of curvature and confinement effects) 
upon the transport via strongly shearing because there are no available experimental data and 
numerical simulations for the same geometric configuration (annular tube with wavy corrugations 
in transverse direction). With a series of forcings : <& = — i? 2 (<%?)/(2To), we can determine the 
enhanced shear rates (dj/dt) due to gravity forcings. From equation (6), we have (up to the 
first order) 

!tL = ^[sinh$ + esin(A;0)-|-cosh$]. (9) 
dt dt R2 

The parameters are fixed below (the orientation effect : sm(k9) is fixed here). r 2 (the mean 

outer radius) is selected as the same as the slip length L Q S . The amplitude of wavy roughness 

can be tuned easily. The effect of wavy-roughness is significant once the forcing (<3?) is rather 

large (the maximum is of the order of magnitude of e[$tanh($)/i?2])- 

If we select a (fixed) temperature, then from the expression of To, we can obtain the shear stress 
r corresponding to above gravity forcings (<£) : 

r = r sinh [sinh(<£) + e sin(fc0) — cosh(<E>)] . (10) 

R2 

There is no doubt that the orientation effect (9) is also present for the amorphous matter. For 
illustration below, we only consider the maximum case : \sm(k6)\ = 1. We shall demonstrate 
our transport results below. The wave number of roughness in transverse direction is fixed to 
be 10 (presumed to be the same for both interfaces of the annular tube) here. 
As the primary interest of present study is related to the possible frictionless transport or for- 
mation of superfluidity (presumed to be relevant to the CGC as mentioned in Introduction) due 
to strong shearing, we shall present our main results in the following. Note that, based on the 
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absolute-reaction-rate Eyring model (of stress-biased thermal activation), structural rearrange- 
ment is associated with a single energy barrier (height) AE that is lowered or raised linearly by 
a (shear) yield stress r. If the transition rate is proportional to the plastic (shear) strain rate 
(with a constant ratio : Co; 7 = Co-Rt, Rt is the transition rate in the direction aided by stress), 
we have 

where uq is an attempt frequency or transition rate, CqUq ~ 70 exp(AE/kBT), or 

It is possible that the frictional resistance (or shear stress) can be almost zero (existence of 
r ~ 0) from above equations (say, equation (11) considering a sudden jump of the resistance). 
The nonlinear character only manifests itself when the magnitude of the applied stress times the 
activation volume becomes comparable or greater in magnitude than the thermal vibrational 
energy. 

Normally, the value of Vh is associated with a typical volume required for a microscopic shear 
rearrangement. Thus, the nonzero transport rate (of the condensed composite (quarks and glu- 
ons) system) as forcing is absent could also be related to a barrier-overcoming or tunneling for 
shear-thinning matter along the wavy-roughness (geometric valley and peak served as poten- 
tial surfaces) in annular tubes when the wavy-roughness is present. Once the geometry-tuned 
potentials (energy) overcome this barrier, then the tunneling (spontaneous transport) inside 
wavy-rough annular tubes occurs. Now, we start to examine the temperature effect. We fix the 
forcing $ to be 1 as its effect is of the order 0(1) for the shear rate. As the gravity forcing 
(SQ) might depend on the temperature (\SQ\ = 2tq&/R2, to = T o(T), V*(= Vh) is presumed 
to be temperature independent here for simplicity). Note that, according to [7], V* = 3V5j/2 
for certain matter during an activation event [7], where V is the deformation volume, ^7 is the 
increment of shear strain. 

As the primary interest of present study is related to the possible phase transition [14-16] or 
formation of superfluidity (presumed to be relevant to the formation of dark matter mentioned 
in Introduction) due to strong shearing, we shall present our main results in the following. We 
performed intensive calculations or manipulations of related physical and geometric parameters, 
considering a hot big-bang universe [14] and examine what happens as it expands and cools 
through the transition temperature T c . The selected temperature range and the activation an- 
ergy follows this reasoning. Note that in unified models of weak and electromagnetic interactions 
T c is of the order of the square root of the Fermi coupling constant [14], G 1 / 2 , i.e. a few hundred 
GeV. Thus the transition occurs when the universe is aged between 1CT 10 and 1CT 12 seconds 
and far above nuclear densities [23]. One possible high-temperature superfluidity formation (or 
nearly frictionless transport) regime [2,4] is demonstrated in Fig. 1. The activation energy 



9 



(AE) is 6 x 10 -10 Joule. In fact, all the results shown in this figure depend on 70 and are thus 
very sensitive to AE (and V/j). Here = 2 and the sudden jump of the shear stress (directly 
linked to the friction) occurring around T ~ 10 12 °K could be the transition temperature for 
the selected AE and Ck- There is a sudden friction drop around two orders of magnitude below 
T ~ 2 x 10 12 °K (V h ~ 4 x KT 11 m 3 ) and it is almost frictionless below T ~ 10 12 °K. The 
temperature regime we identify here is close to the critical (at the critical end point) as well as 
crossover (to the quark-gluon plasma (QGP) at zero baryon chemical potential,) temperature 
mentioned in [2] (cf. [1]). 

The possible reasoning for this formation can be illustrated in Fig. 2. It could be due to the 
strong shearing driven by larger forcings along a confined tube. The shear-thinning (the viscosity 
diminishes with increasing shear rate) reduces the viscosity significantly. One possible outcome 
for almost vanishing viscosity is the nearly frictionless transport. 

Based on the knowledge gained at RHIC, it should be possible to predict experimental results at 
the Large Hadron Collider (LHC), which will collide lead ions at much higher energy densities 
[24]. At this point it's important to predict something about the LHC experiments [13,24] even 
the LHC will bring unanticipated discovery. Fig. 3 illustrates the possible (almost) perfect fluid 
formation or frictionless transport for ultra-high temperature transport of amorphous matter 
(possible color-glass-condensate [24]) considering Pb+Pb collisions at LHC. The activation en- 
ergy is 10" 9 Joule. We can observe there is a sudden friction (shear stress) drop around two 
orders of magnitude below T ~ 7 x 10 17 °K (Vh ~ 2.5 x 10~ 6 m 3 ) and it is almost frictionless 
below T~3x 10 17 °K. 

4 Conclusions 

To conclude in brief, we have obtained critical parameters for the possible perfect fluid formation 
or almost frictionless transport of ultra-high temperature amorphous matter (possible CGC or 
string-like composite (condensed) particles associated with voids) relevant to Au+Au collisions 
at RHIC as well as Pb+Pb collisions at LHC. These critical parameters depend strongly upon 
the temperature, activation energy and activation volume. We shall investigate other relevant 
issues [1,13,25,26] in the future. 
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Figure 1. Comparison of calculated (shear) stresses using an activation energy 6 x 
There is a sharp decrease of shear stress around T~2x 10 12 K. Below 10 12 K, 
the transport of amorphous matter is nearly frictionless (cf. [1,2] for the critical 
temperatures related to Au+Au collisions at RHIC). 
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Figure 2. Increasing shear causes a local energy minimum to flatten until it disappears 
(energy barrier removal or quantum-like tunneling). The structural contribution 
to the shear stress is referred to shear-thinning. 
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Figure 3. Comparison of calculated (shear) stresses using an activation energy 10~ 9 J. 
There is a sharp decrease of shear stress around T ~ 7 x 10 17 K. Below 3 x 10 K, 
the transport of amorphous matter is nearly frictionless (cf. [24] for the possible critical 
temperatures in Pb+Pb collisions at LHC). 



